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Abstract 


Micromechanical  Disk  Array  for  Enhanced  Frequency  Stability  Against  Bias  Voltage 

Fluctuations 
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Professor  Clark  T.-C.  Nguyen,  Chair 

W\g\\-Q  capacitive-gap  transduced  micromechanical  resonators  constructed  via  MEMS 
technology  have  recently  taken  center-stage  among  potential  next  generation  timing  and 
frequency  reference  devices  that  might  satisfy  present  and  future  applications.  Notably, 
oscillators  referenced  to  very  high  Q  capacitive-gap  transduced  MEMS  resonators  have 
already  made  inroads  into  the  low-end  timing  market,  and  research  devices  have  been 
reported  to  satisfy  GSM  phase  noise  requirements  while  only  consuming  less  than  80  pW 
of  power.  Meanwhile,  such  devices  have  also  posted  some  impressively  low  acceleration 
sensitivities,  with  measured  sensitivity  vectors  less  than  0.5  ppb/g. 

Interestingly,  theory  predicts  that  the  acceleration  sensitivity  of  these  devices  should  be 
even  better  than  this,  if  not  for  frequency  instability  due  to  electrical  stiffness.  Indeed, 
electrical  stiffness  is  predicted  to  set  lower  limits  on  not  only  short-term  stability,  but  long¬ 
term  as  well,  especially  when  one  considers  frequency  variations  due  to  charging  or 
temperature-induced  geometric  shifts. 

Pursuant  to  circumventing  electrical  stiffness-based  instability,  this  work  introduces  a 
more  circuit  design-friendly  equivalent  circuit  model  that  uses  negative  capacitance  to 
capture  the  influence  of  electrical  stiffness  on  device  and  circuit  behavior.  This  new  circuit 
model  reveals  that  capacitive -gap  transduced  micromechanical  resonators  can  offer  better 
stability  against  electrical-stiffness-based  frequency  instability  when  used  in  large 
mechanically-coupled  arrays.  Measurements  confirm  that  a  215-MHz  50-resonator  disk 
array  achieves  a  3.5 X  enhancement  in  frequency  stability  against  dc-bias  voltage  variation 
over  a  stand-alone  single  disk  counterpart.  The  new  equivalent  circuit  predicts  the 
measurement  data  and  its  trends  quite  well,  creating  good  confidence  for  using  this  circuit 
to  guide  new  oscillator  and  filter  designs  that,  depending  on  the  application,  can  enhance 
or  suppress  electrical  stiffness. 
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Chapter  1  Introduction 


1.1.  Motivation:  Frequency  Stability  Enhancement  for  Oscillators 

High-0  capacitive-gap  transduced  micromechanical  resonators  constructed  via  MEMS 
technology  have  recently  taken  center-stage  among  potential  next  generation  timing  and 
frequency  reference  devices  that  might  satisfy  present  and  future  applications.  Notably, 
oscillators  referenced  to  very  high  Q  capacitive-gap  transduced  MEMS  resonators  have 
already  made  inroads  into  the  low-end  timing  market,  and  research  devices  have  been 
reported  to  satisfy  GSM  phase  noise  requirements  [1]  [2].  Figure  1.1  presents  a  61 -MHz 
Pierce  oscillator  referenced  to  a  capacitive-gap  transduced  wine  glass  disk  resonator  that 
makes  GSM  phase  noise  requirement  at  both  close  to  carrier  and  far-from-carrier  offsets, 
with  consuming  only  78  microwatts  of  power  [3].  However,  like  any  other  oscillators, 
environmental  fluctuations  such  as  external  vibrations,  bias  voltage  noise,  or  charging,  may 
greatly  degrade  this  performance  [4]  [5]. 

The  work  in  [6]  on  such  capacitive-gap  MEMS  devices  posted  some  impressively  low 
acceleration  sensitivities,  with  measured  sensitivity  vectors  less  than  0.5ppb/g. 
Interestingly,  theory  predicts  that  the  acceleration  sensitivity  of  these  devices  should  be 
even  better  than  this,  if  not  for  frequency  instability  due  to  electrical  stiffness  [6].  Indeed, 


Offset  Frequency  [Hz] 


Figure  1.1.  A  61 -MHz  Pierce  oscillator  referenced  to  a  capacitive-gap  transduced  wine  glass  disk  resonator 
with  its  phase  noise  performance  making  the  GSM  specifications. 
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electrical  stiffness  is  predicted  to  set  lower  limits  on  not  only  short-term  frequency  stability, 
but  long-term  as  well.  As  the  electrical  stiffness  ke  of  a  capacitive-gap  MEMS  resonator  is 
often  determined  by  de-bias  voltage  VP,  electrode-to-resonator  overlap  capacitance  C0,  and 
the  capacitive  gap  spacings  d0 ,  environmental  fluctuations  that  disturb  any  of  these 
parameters  will  cause  instability  in  electrical  stiffness  and  generate  frequency  shift,  as 
shown  in  Figure  1 .2.  For  instance,  theoretical  analysis  of  micromechanical  wine-glass  disk 
resonators  reveals  that  acceleration-induced  changes  in  electrode-to-resonator  gap  spacing 
or  overlap  area  that  in  turn  induce  shifts  in  electrical  stiffness  dominate  among  sources  that 
shift  frequency  during  accelerations  [7].  In  addition,  noise  or  drift  on  the  power  supply 
manifests  as  fluctuations  on  the  resonator  de-bias  Vp  that  obviously  destabilize  the 
electrical  stiffness,  and  thereby,  resonance  frequency  [4]. 

Pursuant  to  suppress  the  electrical  stiffness  and  enhance  frequency  stability,  the  work  in 
this  report  develops  a  new  negative  capacitance  based  equivalent  circuit  model  [8],  which 
reveals  that  capacitive-gap  transduced  micromechanical  resonators  can  offer  better 
frequency  stability  against  environment  fluctuations  when  used  in  large  mechanically- 
coupled  arrays  [9].  The  key  to  enhanced  frequency  stability  is  the  electrode-to-resonator 
capacitance  (Co)  generated  by  the  parallel  combination  of  input/output  electrodes 
overlapping  each  resonator  in  the  array  that  in  turn  reduces  the  efficacy  of  the  bias  voltage- 
controlled  electrical  stiffness.  Here,  the  new  equivalent  circuit  based  on  negative 
capacitance  also  provides  improved  visualization  that  helps  to  identify  methods  to  suppress 
electrical  stiffness  induced  frequency  variation.  The  circuit  model  indicates  that  the  most 
stable  MEMS-based  oscillators  (e.g.,  against  supply  noise  and  acceleration)  are  the 
capacitive-gap  transduced  ones  that  utilize  mechanically-coupled  arrays  of  resonators. 


•  Short-term 


Figure  1.2.  Various  factors  that  can  contribute  to  electrical  stiffness  induced  frequency  instability,  which 
includes  bias  voltage  noise,  environment  vibration,  dielectric  charging,  and  temperature  drift. 
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1.2.  Capacitive-gap  Transduced  Radial  Contour  Mode  Disk  Resonators 


A  capacitive-gap  transduced  micromechanical  resonator  usually  consists  of  a  suspended 
mechanical  structure  and  capacitive  transducers  that  convert  energy  between  electrical 
domain  and  mechanical  domain.  It  generally  offers  the  best  quality  factor  among  different 
types  of  micromechanical  resonators  (which  is  desired  for  low  phase  noise  oscillator 
applications)  [10]  [1 1],  since  they  generally  are  constructed  in  single  high  quality  materials, 
and  thus  suffer  less  from  the  material  interface  losses  that  can  encumber  other  transducer 
types  (e.g.,  piezoelectric).  Figure  1.3  compares  the  quality  factor  of  a  capacitive-gap 
transduced  disk  resonator  and  a  piezoelectric  ring  resonator  at  similar  frequencies  around 
500  MHz  [12]  [13].  The  capacitive-gap  transduced  disk  resonator  achieved  10x  higher  Q 
compared  with  the  piezoelectric  resonator.  Although  piezoelectric  transducers  may  be  more 
successful  in  achieving  lower  impedances  between  50  H  and  377  H  for  matching  to  off- 
chip  wireless  components  [14],  capacitive-gap  transduced  devices  can  scale  down  their 
impedances  by  using  some  of  the  recent  developed  techniques  [15]  [16].  In  addition  to 
better  Q,  capacitive-gap  transduced  resonators  also  offer  more  flexible  geometries  with 
CAD-definable  frequencies,  as  they  can  take  the  form  of  many  shapes,  such  as  beams,  disks, 
rings,  or  plates  [17]  [18]  [10]  [19],  and  their  frequencies  usually  depend  on  lateral 
geometries  that  can  be  defined  by  lithography,  which  enables  the  fabrication  of  multi- 
frequency  devices  on  a  single  chip. 
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Figure  1.3.  Q  Comparison  between  a  capacitive-gap  transduced  disk  resonator  and  a  piezoelectric  ring 
resonator  at  similar  resonance  frequencies. 
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(a)  4r  (b) 


Figure  1.4.  (a)  Schematic  of  a  capacitive-gap  transduced  all-polysilicon  radial  contour  mode  disk  resonator 
in  a  two-port  excitation  and  sensing  configuration;  (b)  Finite  element  simulation  of  radial  contour  mode 
shape. 

Figure  1.4(a)  depicts  a  capacitive  gap  all-polysilicon  contour  mode  disk  type  resonator 
that  will  be  discussed  in  this  work  in  a  typical  bias,  excitation,  and  detection  scheme.  This 
device  comprises  a  2pm-thick  disk  surrounded  by  two  closely  spaced  {do  =  50  nm) 
electrodes  and  supported  by  an  anchored  stem  attached  at  the  center  of  the  disk,  where  the 
radial  displacement  is  minimized  in  radial  contour  mode  shape  shown  in  Figure  1.4(b), 
making  this  attachment  location  a  quasi-nodal  point  (“quasi”,  since  the  stem  has  finite 
radius). The  device  is  excited  into  resonance  via  a  combination  of  a  de-bias  voltage  Vp 
applied  to  the  conductive  polysilicon  resonant  structure  and  an  ac  signal  vt  applied  to  the 
input  electrode,  which  together  induce  a  force  at  the  frequency  of  vt  that  drives  the  disk 
into  resonance  vibration  when  the  frequency  of  xy  matches  the  resonance  frequency.  Once 
vibrating,  the  Fp-biased  time  varying  capacitance  between  the  disk  and  its  output  electrode 
generates  an  output  current  detectable  by  measurement  instrumentation.  The  resonance 
frequency  of  a  radial  contour  mode  disk  resonator  can  be  expressed  by: 


1 

oc  — 
R 


(1) 


where  kr,  mm,  R,  E,  and  p  represent  the  equivalent  stiffness,  dynamic  mass,  disk  radius, 
material  Young’s  modulus,  and  density  of  disk  structure,  respectively. 


1.3.  Traditional  Equivalent  Circuit 

To  conveniently  model  and  simulate  the  behavior  of  a  micromechanical  radial  contour 
mode  disk  resonator  in  a  way  that  facilitates  circuit  analysis,  Figure  1.5  presents  a 
traditional  ac  small-signal  electrical  equivalent  circuit  that  includes  a  core  ler  tank  to  model 
the  mechanical  vibration,  a  variable  capacitor  to  model  the  dependence  of  resonance 
frequency  on  electrical  stiffness  ke  ,  and  a  pair  of  transformers  to  represent  the 
electromechanical  couplings. 
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Figure  1.5.  Classic  equivalent  electrical  circuit  for  a  capacitive  gap  micromechanical  disk  resonator  with 
electrical  stiffness  lumped  into  the  variable  capacitance  cr. 


1.3.1.  Core  ler  Tank 


Despite  its  mechanical  nature,  the  disk  resonator  of  Figure  1.4(a)  still  looks  like  an 
electrical  device  when  looking  into  its  ports,  and  so  can  be  modeled  by  the  electrical 
inductor-capacitor-resistor  ( ler )  equivalents  shown  in  Figure  1.5.  Expressions  for  the 
element  values  in  the  ler  equivalents  take  the  form: 


Ly  cirri’ 


(2) 


where  mm,  kr,  km,  ke  and  bm  are  the  actual  values  of  dynamic  mass,  equivalent  stiffness, 
mechanical  stiffness,  electrical  stiffness,  and  damping  constant  of  the  resonator  being 
modeled  at  the  core  ler  location.  For  radial  contour  mode  shape  with  all  points  on  the 
sidewalls  of  the  disk  moving  the  same  amount,  the  perimeter  of  the  disk  structure  usually 
becomes  the  core  ler  location.  The  equivalent  mass  can  be  obtained  by  dividing  the  total 
kinetic  energy  KEtot  by  one-half  the  square  of  the  velocity  at  the  core  ler  location  [18]. 
Doing  so  on  the  perimeter  of  the  disk  yields: 


KEtot 

T7lm  * 

^ v2(R,9 ) 


2npt  rJl  (hr) dr 

Jl(hK) 


(3) 


where  p,  t,  and  R  are  the  material  density,  thickness,  and  radius  of  the  disk  structure; 
Jn(hr)  is  the  Bessel  function  of  the  first  kind  of  order  n;  and  h  is  a  constant  defined  by 
Young’s  Modulus  E,  density  p,  Poisson  ratio  v  of  the  disk  material,  and  radian  resonance 
frequency  u>0 : 


h  = 


2  E 


2  +  2v 


+ 


Ev 


1  —  v2 


(4) 
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From  (2)  and  (3),  expressions  for  lx,  cr,  and  rx  in  the  ler  tank  at  a  location  on  the  disk 
perimeter  can  be  obtained  using  the  following  relations  [18]: 


lx  =  m 


m 


(5) 


Cr  2 


CL)nmr 


rv  = 


Q 


(6) 

(7) 


1.3.2.  Electromechanical  Coupling  Factor  tjei 


As  the  ler  tank  only  captures  the  mechanical  behavior  of  the  disk  structure,  this 
equivalent  circuit  also  needs  a  transformer  to  model  the  transduction  from  mechanical 
domain  to  electrical  domain  or  vice  versa.  A  succinct  derivation  of  the  transformer  turn 
ratio  rjei,  which  is  also  called  the  electromechanical  coupling  factor,  follows  directly  from 
consideration  of  the  mechanical  forces  generated  by  voltages  applied  to  terminals  of  the 
disk  resonator  depicted  in  Figure  1 .4(a).  A  de-bias  voltage  VP  on  the  resonator  body  and  an 
ac  voltage  signal  vL  at  the  input  together  generate  an  electrostatic  input  force  Ft  in  a  radial 
direction  that  takes  the  following  form: 


1  dCt 

Pi  =  2-W(y'-V'y 


(8) 


where  Ct  represents  the  electrode-to-resonator  overlap  capacitance  when  the  disk  moves  r 
in  radial  direction,  and  dCJdr  is  the  change  in  electrode-to-resonator  overlap  capacitance 
per  unit  radial  displacement  at  each  corresponding  port.  Retaining  only  the  dominant  term 
in  (8)  at  resonance  frequency  yields: 

(dC;\ 

F‘  “  ~Vp\di)Vi  =  ~netVi  (9 


A  first  order  approximation  for  dCJdr  can  be  obtained  from  the  Taylor  expansion  of  the 
Ct  expression  as  follows  [18]: 


QO)  = 


(10) 


which  leads  to  the  expression  of  electromechanical  coupling  factor  rjei  as: 


V  ei 


VpCgi 

d0 


(11) 
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1.3.3.  Frequency  Pulling  Effect  from  Electrical  Stiffness  ke 


Because  the  electrode-to-resonator  capacitance  is  a  nonlinear  function  of  the  disk  radial 
displacement,  as  indicated  in  (10),  there  are  actually  many  more  force  components 
generated  than  represented  in  (9).  In  particular,  Taylor  expanding  (10)  further  to  the  second 
term  yields  the  electrostatic  input  force  Ft  as: 


2VpVi  +  vf  + 


2Vh 


dr 


AVpViV  2vfr N 


dr 


dr 


\C_oi_ 

2  d0 


2Vvvi  + 


2Vp 


dr 


r  2vfr 
-  +  -7^  + 

Clrt 


(12) 


where  the  last  form  includes  only  terms  that  can  generate  components  at  resonance 
frequency.  Inserting  vt  =  FjCOS6L>0t  and  r  =  9fsinu)0t  (where  the  fact  that  the 
displacement  r  is  90°  phase-shifted  from  v;  has  been  accounted  for)  into  (12)  yields: 


Pi 


V  C  ■ 

vp  UOl 

d0 


FjCosu)0t  — 


C 

9fsinu)0t 

d% 


(13) 


At  resonance,  the  second  force  term  in  (13)  is  in  phase  with  the  radial  displacement.  This, 
together  with  the  fact  that  it  is  also  proportional  to  the  displacement  r,  identifies  this  force 
component  as  equivalent  to  a  stiffness,  but  in  this  case,  one  generated  via  electrical  means. 
In  particular,  this  force  component  arises  from  the  increase  and  decrease  in  electric  field 
strength  across  the  electrode-to-resonator  gap  as  the  gap  shrinks  and  grows,  respectively, 
during  mechanical  resonance  vibration.  When  the  disk  sidewall  gets  close  to  the  electrode, 
the  force  pulling  it  into  the  electrode  grows,  and  vice  versa  for  the  other  direction,  i.e.  it 
shrinks  as  the  disk  sidewall  moves  away  from  the  electrode.  Thus,  rather  than  acting  to 
oppose  displacement,  as  is  the  case  for  mechanical  stiffness,  this  force  acts  to  enhance  it, 
which  effectively  makes  it  a  negative  electrical  stiffness  at  port  i,  with  spring  constant  equal 
to: 


k 


ei 


c  ■  c 

uO i  t/2  u 01 

' d i*  pdZ 


(14) 


where  the  last  form  assumes  that  is  much  smaller  than  VP — a  condition  for  which  the 
reader  is  cautioned  isn’t  always  the  case,  such  as  in  micromechanical  mixlers  [20]. 

The  electrical  stiffnesses  of  all  electrodes  will  subtract  from  the  mechanical  spring 
constant  of  the  resonator  at  the  core  ler  location,  changing  the  resonance  frequency  of  a 
two  port  radial  contour  mode  disk  resonator  to: 


Where  kr  is  the  effective  stiffness  of  the  disk  at  any  point  on  its  perimeter,  ke  is  the  total 
electrical  stiffness  contributed  by  all  electrodes  that  surround  the  disk,  km  is  the  purely 
mechanical  stiffness,  and  fnom  is  the  resonance  frequency  of  the  disk  structure  with  zero 
de-bias  voltage  applied.  The  last  expression  in  (15)  uses  the  binomial  expansion  to 
approximate  f0  for  the  case  in  which  the  mechanical  stiffness  is  many  times  larger  than  any 
of  the  electrical  stiffness,  which  is  generally  true  for  devices  in  this  work.  Rearrange  of  (15) 
yields  the  fractional  frequency  change  due  to  electrical  stiffness  as: 


A/_ 

F 

nom 


1  ke 

2  km 


(16) 


1.4.  Outline 

This  report  first  begins  in  Chapter  2  with  a  complete  development  of  the  new  negative 
capacitance  based  equivalent  circuit  to  capture  the  ke  dependent  resonance  frequency  for 
capacitive-gap  radial  contour  mode  disk  resonators,  which  provides  more  insights  into 
methods  that  can  suppress  such  electrical  stiffness  induced  frequency  pulling  effect.  After 
introducing  the  method  of  building  mechanically  coupled  resonator  array  composite  to 
negate  electrical  stiffness  and  enhance  frequency  stability  in  Chapter  3,  Chapter  4 
demonstrates  the  frequency  stability  enhancement  effect  from  micromechanical  disk  array 
composites  with  different  number  of  resonators  and  confirms  the  new  model  accuracy  by 
comparison  with  measured  data  that  includes  plots  of  frequency  versus  de-bias  voltages. 
Finally,  Chapter  5  explores  the  impact  that  use  of  this  equivalent  circuit  should  have  on 
practical  applications. 
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Chapter  2  Negative  Capacitance 
Based  Equivalent  Circuit  Model 


2.1.  Limitations  of  Traditional  Equivalent  Circuit 


Indicated  by  (16),  the  resonance  frequency  of  a  radial  contour  mode  disk  resonator  will 
be  slightly  lower  than  the  mechanical  natural  frequency  of  the  disk  structure,  due  to 
frequency  pulling  effect  of  the  force  generated  by  time-varying  changes  in  electric  field 
strength  as  vibration  changes  the  parallel-plate  capacitive  electrode-to-resonator  gap.  The 
ratio  of  electrical  stiffness  ke  to  purely  mechanical  stiffness  km  actually  captures  the 
magnitude  of  such  frequency  shift.  By  inserting  the  expression  of  parallel-plate  electrode- 
to-resonator  overlap  capacitance  into  (14),  the  total  electrical  stiffness  ke  for  a  two-port 
radial  contour  mode  disk  resonator  can  be  calculated  as: 


where  e  is  the  permittivity  of  the  gap  material  (i.e.,  vacuum  in  this  case)  and  A  is  the  total 
overlap  area  between  the  resonator  and  its  electrodes.  Changes  of  any  variables  in  ke,  such 
as  dc-bias  voltage  noise,  capacitance  variation  due  to  mechanical  vibration,  or  charging 
induced  bias  voltage  drift,  can  cause  frequency  stability  issues  for  capacitive-gap  MEMS 
resonators  [4]  [5]. 

The  traditional  equivalent  circuit  in  Figure  1.5  models  the  influence  of  electrical 
stiffness  on  resonance  frequency  via  the  arrow  through  capacitor  cr  (that  indicates  this 
capacitor  is  tunable)  and  by  setting  the  value  of  cr  equal  to  1  /(km  —  ke).  Although  this 
method  for  capturing  electrical  stiffness  adequately  predicts  the  resonance  frequency,  it 
does  not  convey  clearly  to  a  circuit  designer  the  impact  of  electrical  stiffness  on  the  overall 
circuit  performance.  Modeling  the  electrical  stiffness  in  this  way  in  fact  hides  some  very 
important  capacitive-gap  resonator  behaviors  when  emplaced  into  certain  circuits.  This 
model  also  encourages  designers  to  dismiss  the  impact  of  electrical  stiffness,  since  many 
designers  just  neglect  the  ke  part  in  the  value  of  cr  when  drawing  up  equivalent  circuits. 

To  remedy  the  above  deficiencies  of  the  traditional  equivalent  circuit  model,  this  chapter 
introduces  a  more  circuit  design-friendly  model  presented  in  Figure  2.1  that  captures  the 
influence  of  electrical  stiffness  on  device  and  circuit  behavior  using  a  negative  capacitance 
exactly  equal  in  magnitude  to  the  shunt  static  electrode-to-resonator  overlap  capacitance 
Coi  at  each  electrode  terminal. 
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Port  1  Core  Icr  Port  2 


Figure  2.1.  Negative  capacitance  equivalent  circuit  of  a  two-port  micromechanical  radial  contour  mode  disk 
resonator.  Here,  the  new  equivalent  circuit  uses  a  negative  capacitance  exactly  equal  in  magnitude  to  the 
shunt  static  electrode-to-resonator  overlap  capacitance  C0i  at  each  electrode  terminal  to  capture  the  electrical 
stiffness  induced  frequency  shift. 


2.2.  Negative  Capacitance  Equivalent  Circuit 


To  capture  the  influence  of  electrical  stiffness  more  clearly,  the  new  equivalent  circuit 
needs  to  model  the  electrical  stiffness  and  mechanical  stiffness  separately,  instead  of 
lumping  them  all  together  as  shown  in  Figure  1.5.  In  general,  a  spring  with  stiffness  k  in 
mechanical  domain  can  be  modeled  as  a  capacitor  in  electrical  domain  with  values  of  1/k. 
Parallel  combination  of  stiffness  corresponds  to  series  connection  of  capacitors  in  an 
equivalent  electrical  circuit.  To  attain  the  circuit  of  Figure  2.1,  we  start  with  the  circuit  of 
Figure  1.5  and  first  separate  cr  into  three  capacitors:  cx  =  l/kmto  model  the  mechanical 
stiffness  of  the  resonator  structure,  and  two  negative  capacitors  cei  to  model  electrical 
stiffnesses  generated  at  the  i  port  as  shown  in  Figure  2.2(a),  which  takes  the  following  form: 


c 


ei  — 


d20 

v?coi 


(18) 


where  i  denotes  the  corresponding  port.  Plugging  the  expression  of  micromechanical 
coupling  factor  r\ei  from  (11)  into  (18)  further  yields  the  cei  expression  as: 


d  f 


coi  = 


VpCc 


r  . 

l 

V2ei 


(19) 


Here,  cei  entirely  captures  the  electrical  stiffness,  allowing  cx  in  Figure  2.2(a)  to  represent 
purely  mechanical  stiffness.  Further  reflecting  the  negative  capacitors  cei  through  the 
transformers  on  both  sides  to  outside  the  core  Icr  loop  yields  Cei  in  the  circuit  of  Figure  2.2 
(b),  where  the  physical  shunt  electrode-to-resonator  capacitors  Coi  are  now  matched  by 
series  negative  capacitors  of  exactly  the  same  values,  as  indicated  by  the  following 
expression: 
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Figure  2.2.  Negative  capacitance  small-signal  equivalent  circuits  for  a  two-port  capacitive-gap 
micromechanical  contour  mode  disk  resonator:  (a)  Negative  capacitance  equivalent  circuit  with  electrical 
stiffness  separated  from  mechanical  stiffness,  (b)  Negative  capacitance  equivalent  circuit  with  electrical 
stiffness  reflected  through  transformers  to  outside  the  core  Icr  loop. 

Cei  =  *leicei  =  ~^oi  (20) 

Note  that  the  resonance  frequency  dependency  on  de-bias  voltage  Vp  of  a  radial  contour 
mode  disk  resonator  is  now  entirely  captured  by  the  electromechanical  coupling 
coefficients  rjei.  Thus,  no  arrow  is  needed  through  the  motional  capacitor,  as  needed  in  the 
conventional  circuit  of  Figure  1.5.  More  importantly,  the  value  of  the  core  cx  element  stays 
constant,  as  it  should.  The  negative  Coi  is  also  a  static  capacitor,  just  like  the  positive  Coi  of 
the  physical  shunt  electrode-to-resonator  capacitance.  With  this  circuit,  by  mere  inspection, 
a  designer  can  now  immediately  see  that  the  shunt  Coi  presents  the  opportunity  to 
effectively  negate  the  electrical  stiffness,  suppress  frequency  changes  induced  by  electrical 
stiffness,  and  thereby  stabilize  the  frequency  against  all  variables  in  (17) — something 
highly  desirable  in  some  oscillator  applications. 

2.3.  Mitigation  of  Electrical  Stiffness 

Pursuant  to  determining  which  circuit  configurations  best  promote  ke  cancellation, 
Figure  2.2(b)  adds  load  impedances  ZLi  at  each  terminal.  In  order  for  the  static  electrode- 
to-resonator  capacitor  Coi  to  negate  the  electrical  stiffness  represented  by  —  Coi  , 
approximately  same  amount  of  current  should  flow  through  Coi  and  —Coi,  which  indicates 
that  Coi  should  pass  most  of  the  current  flowing  into  the  parallel  combination  of  Coi  and 
ZLi.  In  other  words,  the  impedance  of  ZLi  should  be  much  larger  than  that  of  Coi,  or: 
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|Zd»— 7“  (21) 

^o^oi 

From  (21),  to  reduce  electrical  stiffness,  both  ZLt  and  Coi  should  be  large.  If  ZLi  has  no 
reactive  component,  operation  at  high  frequency  also  suppresses  electrical  stiffness. 
However,  increasing  operating  frequency  may  conflicts  with  application  requirement 
because  the  resonance  frequency  is  usually  application  oriented,  which  means  a  designer 
may  not  have  the  freedom  to  change  it.  In  addition,  as  a  resonator  might  be  used  with 
reconfigurable  drive  and  sense  circuits  that  present  non-constant  loads,  it  would  be  also 
challenging  to  achieve  better  frequency  stability  by  purely  relying  on  optimizing  load 
impedances.  These  two,  together,  lead  to  the  only  practical  choice  to  enhance  frequency 
stability  against  electrical  stiffness,  which  is  to  boost  a  resonator’s  electrode -to-resonator’s 
overlap  capacitance  C0,  e.g.,  utilizing  resonators  with  solid  dielectric  gaps  [21]  [22],  or 
building  array  of  resonators. 

The  work  in  this  report  focuses  on  coupled-array-based  methods  to  increase  Coi  towards 
better  frequency  stability  against  changes  in  electrical  stiffness.  In  addition,  the  arraying 
strategy  also  reduces  series  motional  resistance  without  degrading  linearity.  This  is  in 
contrast  to  other  available  approaches  to  increase  electrode -to-resonator  capacitance  Coi, 
such  as  reducing  the  electrode-to-resonator  gap  spacing  ,and  raising  the  effective 
permittivity,  all  of  which  eventually  begin  to  degrade  the  device’s  linearity. 
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Chapter  3  Micromechanical  Disk 
Array  Composite 


3.1.  Schematic  of  a  Radial  Contour  Mode  Disk  Array  Composite 


Figure  3.1(a)  presents  a  micromechanical  disk  array-composite  constructed  by 
mechanically  linking  individual  disk  resonators  via  half-wavelength  coupling  beams.  Here, 
coupling  of  resonators  yields  a  coupled  multi-mode  system,  where  at  each  mode,  all 
resonators  vibrate  at  precisely  the  same  frequency  [23],  allowing  their  outputs  to  be 
combined  to  boost  input  and  output  currents,  thereby  decreasing  the  motional  resistance 
and  increasing  power  handling.  As  first  demonstrated  in  [1],  the  use  of  half-wavelength 
couplers  both  insures  that  all  resonators  vibrate  in  phase,  as  shown  in  Figure  3.1(b);  and 
spreads  the  mode  frequencies  apart,  making  it  easier  to  select  a  specific  mode  (when  only 
one  is  wanted,  e.g.,  in  an  oscillator  application)  by  proper  electrode  phasing. 

As  described  in  [16],  the  mechanical  connection  of  resonators  in  Figure  3.1(a)  actually 
realizes  a  multi-pole  filter  structure  that  now  has  several  modes  of  vibration.  Each  modal 
peak  corresponds  to  a  state  where  all  resonators  are  vibrating  at  exactly  the  same  frequency. 
Figure  3.2  presents  finite  element  simulations  showing  the  different  modes  of  similar 
structure  with  3  constituent  resonators,  which  are  distinguishable  by  the  phasings  between 
the  resonators. 


Figure  3.1.  (a)  Schematic  of  a  radial  contour  mode  disk  array-composite  with  disks  linked  by  X/2  coupling 
beams  to  enforce  in  phase  vibration  of  each  individual  resonator,  (b)  Simulated  modal  depiction  of  a  3- 
resonator  disk  array  with  all  disks  vibrating  in  phase  in  radial-contour  mode  shapes. 
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To  only  select  the  desired  in-phase  mode  for  a  micromechanical  disk  array,  it  is 
advantageous  to  first  separate  other  unwanted  modes  as  far  apart  as  possible.  Since  the 
bandwidth  of  a  mechanical  filter  is  proportional  to  the  stiffness  of  its  resonator-to-resonator 
coupling  springs,  the  first  step  in  selecting  a  single  mode,  while  suppressing  others,  is  to 
couple  the  resonators  with  very  stiff  springs.  Ideally,  coupling  beams  with  lengths 
corresponding  to  half-wavelength  of  acoustic  waves  propagating  in  the  resonator  will  shift 
the  frequency  of  unwanted  modes  to  infinity. 

In  general,  a  coupling  beam  can  be  modeled  as  an  acoustic  transmission  line — the 
mechanical  analog  to  the  familiar  electrical  transmission  line  that  takes  the  following 
expression: 


F  i 

'cos(al)  jZmsin(aiy 

r  p 

ri 

u 

— 

jsinial ) 

— — -  costal) 

r2 

u 

where  Ft  and  Xi  are  the  force  and  velocity  at  corresponding  ports,  l  is  the  length  of  the 
coupling  beam.  Zm  and  a  in  (22)  are  acoustic  characteristic  impedance  and  acoustic  wave 
propagation  constant,  respectively,  which  take  the  forms: 


cc  — — ,  Zm  Asyj  Ep 

•JeJp 


(23) 


where  u>  is  the  radian  frequency  of  the  acoustic  wave,  p  is  the  material  density,  E  is  the 
Young’  modulus  of  the  beam  material,  and  As  is  the  beam  cross  sectional  area. 

Figure  3.3  shows  the  duality  between  an  acoustic  coupling  beam  and  an  electrical 
transmission  line.  Because  the  two  port  modeling  of  an  electrical  transmission  line  can  be 
simplified  by  using  an  impedance  T-network,  a  similar  transformation  can  also  apply  on 
the  mechanical  acoustic  transmission  line,  which  generates  a  simplified  impedance  T- 
network  model  as  shown  in  Figure  3.4  [23],  with  series  impedance  za  and  shunt  impedance 
zc  taking  the  following  forms: 


Electrode 


Figure  3.2.  FEA  simulated  mode  shapes  of  a  mechanically  coupled  three -resonator  radial  contour  mode  disk 
array. 
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Figure  3.3.  Two-port  modeling  of  a  mechanical  acoustic  coupling  beam  in  analogy  to  an  electrical 
transmission  line. 


Figure  3.4.  T-network  modeling  of  a  mechanical  acoustic  coupling  beam  with  the  force  (F,)  in  analogy  to  a 
voltage  in  electrical  domain  and  the  velocity  (X,)  in  analogy  to  a  current  in  electrical  domain. 


za  =  jAsjEp  tan  ( 


0) 


\JeTp 


(24) 


(25) 


The  beam  length  of  a  half- wavelength  coupling  can  be  expressed  as: 


l 


s 


\  _JeTp  _  4Wp 

2  2f  71  a) 


(26) 


Plugging  (26)  into  (24)-(25)  yields  \zc\  =  oo  and  \za\  =  0,  which  makes  the  shunt  path 
of  Figure  3.4  an  equivalent  open  and  forces  the  same  velocity  at  the  two  ends  of  the  beam. 
It  means  that  all  resonators  can  only  vibrate  in  phase  with  same  amplitude,  which 
effectively  knocks  out  all  other  unwanted  modes.  It  should  be  noted  that  the  use  of  half 
wavelength  coupling  beams  serves  to  spread  the  modal  peaks  of  the  filter  structure  in 
Figure  3.2  away  to  infinity,  which  facilitates  the  selection  of  one,  and  only  one,  of  the 
modes. 

In  general,  fabrication  variations  and  design  errors  exist,  which  may  not  be  able  to 
generate  exact  half- wavelength  coupling  beams  and  thereby  can  only  separate  the  peaks  by 
finite  distance  in  frequency  domain.  However,  because  each  mode  exhibits  a  unique 
resonator  phasing,  a  single  mode  can  still  be  selected  by  choosing  the  input  ac  signal  and 
output  electrode  configuration  to  match  the  phasing  of  the  desired  mode,  e.g.  directly 
connecting  all  output  electrodes  together  will  effectively  cancel  the  output  current  from  the 
mode  shape  in  which  the  center  resonator  keeps  still  and  the  two  side  resonators  vibrate 
out-of-phase,  as  shown  in  Figure  3.2. 
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Once  a  single  mode  is  selected,  the  structure  practically  behaves  as  a  single  resonator, 
but  with  a  current  handling  ability  equal  to  the  sum  of  the  currents  from  all  constituent 
resonators.  Thus,  an  TV-resonator  array  can  handle  N times  more  power  and  achieve  N times 
smaller  motional  resistance  than  a  single  resonator. 

3.2.  Frequency  Stability  Enhancement  of  a  Disk  Array  Composite 

Of  course,  this  coupled  array  strategy  not  only  produces  a  larger  total  output  current  and 
smaller  motional  resistance  [16],  but  also  attains  much  larger  electrode -to-resonator 
overlap  capacitors  Coi,  and  thus  reduces  the  electrical  stiffness  and  its  associated  instability, 
according  to  (21). 

Because  a  micromechanical  disk  array-composite  behaves  like  a  single  disk  resonator 
when  only  one  in-phase  mode  is  selected,  its  equivalent  circuit  uses  the  same  schematic 
and  topology  as  that  of  a  single  disk  resonator,  as  shown  in  Figure  3.5(a).  For  a  disk  array 
with  N  mechanically  coupled  disk  resonators,  the  stiffness  km,  mass  mm,  damping  bm, 
electromechanical  coupling  factor  rjei,  and  electrode-to-resonator  overlap  capacitance  Coi 
are  all  N  times  larger  than  for  a  single  disk  resonator.  Thus,  simple  multiplication  or 
division  by  N  is  all  that  is  needed  to  derive  coupled  array  equivalent  circuit  element  values 
from  those  of  a  single  disk  resonator,  as  shown  in  Figure  3.5(b).  Note  that  although  the 
electrical  stiffness  goes  up  by  N  according  to  (17),  so  does  the  mechanical  stiffness,  so  their 
ratio  ke/km  remains  the  same  and  the  fundamental  efficacy  by  which  ke  pulls  the  frequency, 
as  governed  by  (15),  does  not  increase.  Instead,  with  Coi  presenting  a  much  smaller 
impedance,  more  current  flows  through  Coi  than  Zlu  allowing  it  to  cancel  more  of  Cei, 
thereby  negating  the  electrical  stiffness  induced  frequency  instability  via  circuit  interaction. 

(a)  Port  1  Core  Icr  Port  2 

A  _ _ A _ _  A 


(b) 


Port  1  Core  Icr  Port  2 

A  _ _ A _ _  A 


Figure  3.5.  (a)  Negative  capacitance  equivalent  circuit  of  a  disk  array-composite.  With  only  the  in-phase 
mode  selected,  the  structure  practically  behaves  like  a  single  resonator,  making  it  share  the  same  equivalent 
circuit  topology  as  a  single  disk  resonator,  (b)  Negative  capacitance  equivalent  circuit  of  a  disk  array  with  N 
resonators  based  on  element  values  in  the  single  resonator  equivalent  circuit  in  Figure  2.2(b) 
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3.3.  Design  Example 

Table  I  summarizes  the  needed  expressions  (with  reference  to  Figure  3.5(a)  for  variables) 
while  also  succinctly  presenting  a  design  flow  to  achieve  a  radial-contour  mode  disk  array 
with  N  resonators  like  that  of  Figure  3.1(a)  with  a  specific  operating  frequency  and  with 
specific  motional  resistance  Rx  between  its  electrode  terminals.  The  work  in  this  report 
also  designs  several  215-MHz  (R  =  12.8  pm)  radial  contour  mode  disk  arrays  with 
various  number  of  mechanically  coupled  resonators  (N=  1,8,16,  and  50),  whose  frequency 
stability  enhancement  effect  will  be  analyzed  in  Chapter  4.  Table  II  summarizes  the  design 
parameters,  the  equivalent  circuit  variables,  and  the  electrical  stiffness  induced  frequency 
shift  for  four  different  type  of  disk  array  designs  which  have  the  same  disk  radius  (R  = 
12.8  pm)  and  coupling  beam  length,  but  different  number  of  resonators. 


Table  I:  Radial  contour  mode  disk  array  design  equations  and  procedure  summary 


Objective/Procedure 

Parameter 

Relevant  Design  Equations  for  a  Given  Parameter 

Eq. 

Contour  Mode  Air  Gap  da 

Disk  0  =  0;  2  / 

Electrode  Electrode 

0  =  0;i 

Solve  For 
b 

*  MS)  , 

s*his>  =  1~v 

(27) 

Angular 

Resonance 

Frequency 

^ nom 

8  (  E 

“nom-  r  L(  1_V2) 

(28) 

Given:  co„nm,  VP,  N, 

Rxij  (=  resistance 

between  terminals  i  and 

./)• 

Find:  R ,  dn. 

1.  Choose  E,  p,  and  v 
by  choice  of 
structural  material. 

2.  Choose  thickness  t. 

3.  Use  (28)  to  find  the 

R  needed  to 
achieve  co,wm ■  Use 
(27)  to  get  8  in  the 
process. 

4.  Use  (32)  to  find  the 
d0  needed  to 
achieve  RxlJ. 

5.  (29),  (30)  and  (31) 
yield  all  needed 
values  in  the  ac 
small  signal 
equivalent  circuit 
in  Figure  3.5(a). 

Core 

Equiv. 

Circuit 

Elements 

2npt  j'{  rjl  (hr)  dr 

L„  =  Nmm  =  N - 5— — - , 

xn  m  Jl(hR) 

,  I  ^0 nomP 

(  2E  +  ^  ) 

1^2  +  2v  1  l-v2j 

1  1  ...  ^nomNlTlm 

C™~  Nkm~  r™-N<-  Q 

(29) 

Static 

Overlap 

Capacitor 

n  _  KJ6o(Si2- 0tl)Rt 

^ oni  ™  1 

d0 

(30) 

Electro¬ 

mechanical 

Coupling 

Coefficient 

r  ■ 

ir  u om 

Veni 

d0 

(31) 

Motional 

Resistance 

„  Un 

XlJ  Venihenj 

(32) 
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Table  II:  Comparison  of  design  variables  and  equivalent  circuit  element  values  for  215-MHz  radial 
contour  mode  disk  arrays  with  different  number  of  mechanically-coupled  resonators. 


Parameter 

Source 

Number  of  Resonators 

Units 

Design  Variables 

Number  of  Resonators,  N 

N/A 

1 

8 

16 

50 

N/A 

Disk  Radius,  R 

Layout 

12.8 

12.8 

12.8 

12.8 

pm 

DC  Bias  Voltage,  VP 

N/A 

9 

9 

9 

9 

V 

Disk  Thickness,  t 

Fabrication 

2 

2 

2 

2 

pm 

Electrode-to-Resonator 
Gap,  d0 

Fabrication 

50 

50 

50 

50 

nm 

Natural  Frequency,  fnom 

Equation  (28) 

215 

215 

215 

215 

MHz 

Coupling  Beam  Length  l 

Layout 

21.2 

21.2 

21.2 

21.2 

pm 

Electrode  Span  Angle,  8 

Layout 

133 

133 

133 

133 

(°) 

Quality  Factor,  Q 

Measurement/ 

Estimate 

10,000 

10,000 

10,000 

10,000 

N/A 

Equivalent  Circuit 
Elements 

Equivalent  Inductance, 

Ln 

Equation  (29) 

1.81xl0"12 

1.45X1011 

2.89xl0‘n 

9.03xl0n 

H 

Equivalent  Capacitance, 

cxn 

Equation  (29) 

3.03xl07 

3.79xl0"8 

1.89x10  s 

6.06xl09 

F 

Equivalent  Resistance, 

Ten 

Equation  (29) 

2. 44x1 0'7 

1.95xl0"6 

3.91xl06 

1.22xl0'5 

a 

Overlap  Capacitance, 
Cmi  —  Cn2 

Equation  (30) 

9.2 

73.8 

148 

462 

fF 

Electromechanical 
Coupling,  Tjenl  =  Ven2 

Equation  (31) 

1.46xl0"6 

1.17x10  s 

2.33xl05 

7.29xl0-5 

N/A 

Electrical  Stiffness 
Induced  Frequency  Shift 
(Jo  fnom)  /  fnom 

Equivalent 
Circuit  of 
Figure  3.5(b) 

-71.3 

-56.6 

-38.7 

-19.8 

ppm 

3.4.  Fabrication  Process 

The  fabrication  process  for  the  all-polysilicon  contour  mode  disk  resonator  arrays  of 
this  work  deviates  from  previous  ones,  such  as  that  of  [18],  in  that  it  does  not  use  self- 
aligned  peg-stem  anchors  and  it  employs  chemical  mechanical  polishing  (CMP)  to  remove 
electrode  overhangs.  The  process  begins  with  film  depositions  and  etches  identical  to  those 
of  [18]  to  achieve  the  substrate  isolation  layer,  polysilicon  interconnects,  and  the  bottom 
sacrificial  layer.  At  this  point,  unlike  previous  self-aligned  processes,  a  mask  is  used  to 
define,  pattern,  and  etch  stem  anchor  holes  into  the  bottom  sacrificial  oxide,  followed  by  a 
2-pm  LPCVD  in-situ  doped  polysilicon  film  that  fills  the  holes  to  form  the  stems  and  serves 
as  the  resonator  structural  material.  Here,  an  AMSL300  DUV  Stepper  is  used  to  realize 
very  precise  alignment,  with  less  than  lOOnm  error — good  enough  to  achieve  Q’ s 
comparable  to  those  of  devices  with  self-aligned  stems. 

After  depositing  an  oxide  hard  mask  over  the  structural  polysilicon,  disk  devices  and 
coupling  links  are  patterned  and  etched  as  before,  the  gap-defining  sacrificial  sidewall 
oxide  spacer  is  deposited,  electrode  to  interconnect  contact  vias  are  etched,  and  the 
polysilicon  electrode  material  LPCVD’ ed  3 -pm  thick,  all  to  yield  the  cross-section  of 
Figure  3.6(a).  At  this  point,  the  process  again  deviates  from  that  of  [18]  in  that  before 
patterning  and  etching  the  electrodes,  the  top  poly  silicon  is  first  CMP’ed  down  to  the  hard 
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Figure  3.6.  Cross-sections  showing  the  last  few  steps  in  the  fabrication  process  for  a  215-MHz  all-polysilicon 
disk  resonators  with  CMPed  electrodes  (a.  b),  and  a  comparison  to  previous  work’s  fabrication  process,  (c). 

mask.  This  step  removes  the  electrode  overhangs  of  Figure  3.6(a),  achieving  the  final  cross 
section  of  Figure  3.6(b).  When  compared  to  the  previous  cross  section  of  [11],  cf.  Figure 
3.6(c),  this  new  process  greatly  increases  the  pull-in  voltage  of  these  devices,  where  contact 
between  the  disk  and  the  electrode  overhang  is  often  the  first  to  occur  when  de-bias 
voltages  increase.  As  before,  devices  are  released  in  49  wt.  %  hydrofluoric  acid  for  ~40min. 

Figure  3.7  and  Figure  3.8  present  SEMs  of  fabricated  215-MHz  50-nm  capacitive-gap 
transduced  half- wavelength-coupled  disk  arrays  employing  16  and  50  resonators, 
respectively,  in  straight  line  and  rectangular  placement  configurations. 


Figure  3.7.  SEM  of  a  215-MHZ  50nm  capacitive-gap  contour  mode  disk  array  employing  16  resonators. 
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Figure  3.8.  SEM  of  a  215-MHZ  50nm  capacitive-gap  contour  mode  disk  array  employing  50  resonators. 
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Chapter  4  Measurement  Results 

To  gauge  the  degree  to  which  arraying  enhances  frequency  stability  against  de-bias 
voltage  changes,  measured  plots  of  frequency  versus  de-bias  voltage  are  in  order,  for  both 
mechanically  coupled  disk  array-composite  resonators  using  various  numbers  of  resonators, 
as  well  as  for  a  single  stand-alone  disk  for  comparison.  To  this  end,  measurements  were 
made  under  a  2ptorr  vacuum  environment  in  a  Lakeshore  FWPX  vacuum  probe  station  and 
using  an  Agilent  E5071C  network  analyzer  in  the  direct  two-port  excitation  and  sensing 
scheme.  Figure  4. 1  presents  the  direct  two-port  measurement  setup  used  in  this  work  for  a 
radial  contour  mode  disk  resonator,  which  records  the  S21  parameter  of  the  devices  to 
observe  the  resonance  frequency  as  a  function  of  the  de-bias  voltage  VP. 

4.1.  Frequency  Stability  Against  DC-bias  Voltages 

Figure  4.2  plots  the  measured  S21  parameter  of  a  single  radial  contour  mode  disk 
resonator  under  different  de-bias  voltage  conditions  (Vp  =  5  V,  10  V  and  15  V).  It  observes 
a  180  ppm  frequency  change  as  Vp  varies  from  15  V  down  to  5  V,  which  clearly  shows  the 
electrical  stiffness  induced  frequency  pulling  effect.  It  also  indicates  that  the  frequency 
stability  of  a  single  disk  resonator  will  be  susceptible  to  all  factors  that  changes  electrical 
stiffness,  such  as  environment  acceleration  or  de-bias  voltage  noise. 


Figure  4.1.  S21  direct  measurement  setup  for  a  polysilicon  215-MHz  capacitive-gap  radial-contour  mode 
disk  resonator. 
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Figure  4.2.  Frequency  spectrum  of  a  single  radial  contour  mode  disk  resonator  under  different  dc-bias  voltage 
conditions.  (Vp  =  5  V,  10  V,  and  15  V) 

To  demonstrate  the  frequency  stability  enhancement  against  dc-bias  voltages  from 
micromechanical  disk  array  composite,  Figure  4.3  plots  the  fractional  frequency  change 
along  the  v-axis  against  dc-bias  voltages,  for  radial  contour  mode  disk  arrays  with  various 
number  of  resonators,  e.g.,  N  =  1,  8,  16,  and  50.  The  measured  curves  clearly  show  a 
shrinking  frequency  dependence  on  dc-bias  voltage  as  the  number  of  resonators  used  in  an 
array  increases.  In  particular,  the  50-resonator  215-MHz  disk  array  experiences  a  20  ppm 
frequency  change  when  Vp  varies  over  a  7  V  span,  from  2  V  to  9  V,  3.5 X  smaller  than  the 
70  ppm  of  a  stand-alone  device. 
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Figure  4.3.  Measured  curves  of  resonance  frequency  versus  dc-bias  voltage  Vp  plotted  against  simulation 
using  negative  capacitance  equivalent  circuit  models  for  disk  arrays  with  N=  1,  N=  8,  N=  16,  and  N=  50. 
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Table  III:  Electrode-to-resonator  overlap  capacitance  and  trace  resistances  for  radial  contour  mode  disk 

array  with  different  number  of  resonators 


Resonator  Number  (TV) 

Overlap  Capacitance  Cc,  [fF] 

Ru,tot  LkOJ 

1 

9.2 

1.56 

8 

73.8 

4.04 

16 

148 

5.24 

50 

462 

2.92 

To  confirm  the  validity  of  the  negative  capacitance  equivalent  circuit  of  Figure  3.5(b), 
simulated  plots  using  this  circuit  are  also  included  in  Figure  4.3,  showing  very  good 
agreement  between  theory  and  measurement.  These  simulations  assume  the  load 
impedance  ZLt  derives  from  a  combination  of  series  trace  resistance  RL1  and  RL2  from 
input/output  electrode  leads,  as  well  as  similar  lead  resistance  RL3  from  the  VP  port,  all 
shown  in  Figure  4.1.  Since  probe  coax  and  bond  pad  capacitance  are  nulled  by  calibration, 
the  ZLi  in  these  measurements  is  mainly  resistive.  It  should  be  noted  that  the  total 
equivalent  load  resistance  RLi)tot  °f  a  disk  array  in  this  work  is  generally  larger  than  that 
of  a  single  resonator  and  actually  increases  with  the  number  of  disks  in  the  array.  This 
comes  about  because,  at  least  in  the  current  layout,  the  distance  between  the  Vp  pad  and  the 
furthest  resonator  increases  as  the  number  of  resonators  increases,  as  seen  from  Figure  3.7, 
resulting  in  a  corresponding  increase  in  series  resistance  RL3.  Table  III  summarizes  overlap 
capacitance  and  equivalent  load  resistance  values  for  disk  arrays  with  N=  1,  8,  16,  and  50, 
clearly  showing  larger  values  as  the  number  of  resonators  increases,  at  least  for  straight 
line  arrays  with  7V=8  and  N=  16,  like  that  of  Figure  3.7.  The  series  resistance  of  the  50- 
resonator  array  depicted  in  Table  III  is  actually  smaller  than  those  of  the  8-  and  1 6-resonator 
ones,  since  its  layout  uses  a  rectangular  or  matrix  topology,  rather  than  a  straight  line,  so 
the  average  distance  of  its  resonators  from  its  Vp  pad  ends  up  being  smaller. 

According  to  (21),  increases  in  load  resistance  like  those  in  Table  III  should  also 
contribute  to  an  overall  nulling  of  the  electrical  stiffness,  and  thereby  enhance  frequency 
stability  against  dc-bias  voltage  fluctuations.  Indeed,  as  the  number  of  array  resonators 
increases,  electrical  stiffness  erodes  due  to  increases  in  both  electrode-to-resonator  overlap 
capacitance  and  load  resistance — a  double  whammy  effect  perfectly  predicted  by  the 
negative  capacitance  equivalent  circuit. 


4.2.  Motional  Resistance 

Figure  4.4  presents  the  de-embedded  frequency  characteristics  measured  using  two-port 
configurations  under  vacuum  for  a  single  radial  contour  mode  disk  resonator  and  a  50- 
resonator  disk  array.  To  allow  for  direct  comparison  of  motional  resistance,  the  same  dc- 
bias  voltage  was  applied  to  each  device  for  measurement.  The  50-resonator  disk  array 
achieves  a  motional  resistance  of  2.4  kQ,  which  is  a  28x  improvement  over  the  single  disk 
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Figure  4.4.  Frequency  spectrum  comparison  of  a  50-resonator  radial  contour  mode  disk  array  with  a  single 
disk  resonator,  both  are  measured  with  same  dc-bias  voltage. 

resonator.  The  discrepancy  between  the  measured  28x  improvement  and  the  predicted  50x 
reduction  in  Rx  may  be  attributable  to  the  fact  that  the  50-resonator  disk  array  has  a  lower 
Q  of  10,965,  compared  with  the  single  resonator’s  Q  of  16,664.  Because  of  lower  motional 
resistance,  the  50-resonator  disk  array  will  suffer  more  Q  loading  effect  from  the  trace 
resistance,  which  leads  to  a  lower  quality  factor.  By  taking  into  account  the  Q  difference, 
the  50-resonator  disk  array  actually  achieves  45 x  smaller  motional  resistance  compared 
with  single  disk  resonator,  which  is  very  close  to  the  ideal  value  of  50. 
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Chapter  5  Conclusions 

5.1.  Impact  on  Applications 

The  new  equivalent  circuit  based  on  negative  capacitance  provides  improved 
visualization  that  helps  to  identify  methods  to  suppress  electrical  stiffness  induced 
frequency  variation. 

A.  Reference  Oscillator  Design  Insights 

From  a  practical  circuit  perspective,  the  use  of  small  source  and  load  resistances  equates 
to  the  use  of  voltage  drive  and  current  sensing.  Conversely,  the  use  of  large  source  and  load 
resistances  equates  to  the  use  of  current  drive  and  voltage  sensing.  From  the  discussion  in 
Chapter  3,  voltage  drive  makes  most  sense  where  frequency  tuning  is  needed,  such  as  for 
a  tunable  filter  passband  application,  or  a  voltage-controlled  oscillator.  On  the  other  hand, 
current  drive  is  most  appropriate  when  frequency  stability  is  paramount. 

One  example  application  where  the  choice  of  drive  and  sense  type  could  make  a  big 
difference  is  that  of  a  reference  oscillator  that  must  be  stable  against  a  variety  of 
environmental  perturbations.  These  include  acceleration,  power  supply  noise,  drift,  and 
undue  charging  (e.g.,  due  to  radiation),  all  of  which  can  induce  instability  in  the  oscillator’s 
frequency.  The  most  significant  mechanism  for  frequency  instability  caused  by  these 
particular  perturbations  ends  up  being  instability  in  the  electrical  stiffness.  As  revealed  by 
the  negative  capacitance  equivalent  circuit  in  Figure  2.2(b),  frequency  dependence  on 
electrical  stiffness  can  essentially  be  nulled  by  using  a  high  impedance  input,  high 
impedance  output  sustaining  amplifier.  Here,  high  impedance  is  defined  relative  to  the 
impedance  presented  by  the  shunt  Coi  at  the  resonator  I/O  terminals,  i.e.,  the  resistance 
presented  by  the  drive/sense  circuit  loading  each  I/O  port  is  considered  “high”  when  it  is 
at  least  5  times  larger  than  1  /(sCoi).  This  suggests  that  to  maximize  frequency  stability 
against  environmental  variations,  a  Pierce  oscillator  [24]  configuration  would  be  a  better 
choice  than  the  commonly  used  transresistance  sustaining  amplifier. 

On  the  other  hand,  if  frequency  tuning  is  important,  and  slight  instabilities  due  to 
environmental  perturbations  can  be  tolerated,  a  low  impedance  input,  low  impedance 
output  sustaining  amplifier  is  most  suitable,  such  as  a  transresistance  amplifier  [25]. 

B.  Device  Design  Insights 

For  the  case  where  frequency  stability  is  of  most  interest,  the  negative  capacitance 
equivalent  circuit  offers  a  very  important  insight:  The  highest  stability  against  changes  in 
bias  voltage  or  overlap  capacitance  comes  when  the  positive  Coi s  can  cancel  the  negative 
ones.  It  will  ensure  maximum  stability  if  Coi  can  be  made  much  greater  so  that  its 
impedance  is  smaller  than  any  load  resistance  on  a  given  terminal. 
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This  simple  fact  now  reveals  several  methods  by  which  the  frequency  stability  of  a 
capacitive-gap  transduced  resonator  can  be  maximized: 

1)  Utilize  large  arrays  of  resonators,  perhaps  mechanical  coupled  into  array  composites 
like  that  of  [16]  and  [9].  The  use  of  many  devices  increases  the  total  input  or  output 
Coi  so  that  it  swamps  any  terminal  load  resistance. 

2)  Utilize  solid  dielectric  gaps,  such  as  demonstrated  in  [21]  and  [22].  Here,  simply 
raising  the  gap  permittivity  again  raises  Coi  so  that  it  dominates  over  any  terminal 
resistance. 

3)  Employ  small  gap  spacings,  which  not  only  improves  motional  resistance,  but  also 
raises  Coi  to  cancel  more  of  the  electrical  stiffness  that  is  represented  by  -Coi. 

Needless  to  say,  design  insights  like  these  made  possible  by  the  negative  capacitance 
equivalent  circuit  that  can  greatly  improve  performance  or  enhance  specific  performance 
parameters  over  others,  are  quite  invaluable.  There  are  sure  to  be  many  other  examples, 
from  oscillators,  to  filters,  to  mixlers  [20],  where  explicit  representation  of  the  negative 
capacitance  in  a  micromechanical  resonator’s  equivalent  circuit  proves  instrumental  to 
maximizing  performance. 

5.2.  Conclusion 

By  explicitly  modeling  electrical  stiffness  in  a  capacitive-gap  transduced 
micromechanical  resonator  via  a  negative  capacitance,  the  equivalent  circuit  of  this  work 
provides  designers  with  needed  visual  cues  that  enable  them  to  more  intelligently  choose 
device  topologies  to  maximize  the  performance  of  resonator  applications,  especially  where 
frequency  stability  is  paramount.  The  50-resonator  capacitive-gap  transduced 
micromechanical  disk  array  demonstrated  here  enables  a  mere  20ppm  frequency  change 
over  a  7V  de-bias  voltage  variation — a  3.5 X  reduction  over  a  single  stand-alone  disk 
resonator  counterpart.  That  the  new  negative  capacitance  equivalent  circuit  model 
introduced  perfectly  predicts  this  phenomenon,  while  also  aiding  circuit  visualization, 
bodes  well  for  its  continued  use  in  future  resonator  circuits  for  which  tailored  electrical 
stiffness  strengths  are  desired.  Indeed,  the  demonstrated  stability  enhancing  attributes  of 
mechanically-coupled  arrays  that  make  them  less  vulnerable  to  de-bias  voltage  noise, 
dielectric  charging,  and  external  vibrations,  together  with  already  demonstrated  array- 
derived  reductions  in  the  standard  deviation  of  array  resonance  frequency  [26],  present 
strong  cases  for  a  more  prevalent  use  of  arrays  in  next  generation  MEMS-based  frequency 
reference  devices. 
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